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By integrating a set of model equations for Joseph- 
son ladders subjected to a uniform transverse bias 
current we have found almost all of the kinds 
of breathers described in recent experiments, and 
closely reproduced their voltage-current characteris- 
tics and switching behaviour. Our main result is that 
a chaotic transient occurs in the switching process. 
The growth of tiny perturbations during the chaotic 
transient causes the new breather configuration to be 
extremely sensitive to the precise history of the ini- 
tial breather and can also cause the new breather to 
have a new centre of symmetry. 



The recent discovery |^, |3 of roto-breathers in 
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Josephson ladder arrays (Fig. If) driven by a trans- 
verse dc bias current have shown not only that these 
localised excitations exist but that they exhibit re- 
markable behaviour. 

The Josephson ladder has been studied theoreti- 
cally for many years. In the absence of a driving cur- 
rent the interaction between vortices has been found 
to be exponential and this leads to the vortex density 
exhibiting a devil's staircase as the magnetic field is 
increased jl], ||. Quantum fluctuations ||, meta-stable 
states and inductance effects have been studied 
as has the interaction of vortices with a transverse dc 
bias current H fl0| |. 

A roto-breather|ll| is a stable group of vertical 
junctions rotating together (6j — 6j ~ u>t, see Fig. 0). 
Such solutions have been studied theoretically (l^] in 
ladders driven by a uniform transverse ac bias cur- 
rent and are also stable |l3| for the case of a dc 
bias current. Rotation at a single vertical junction 
was studied numerically [fl4| |l5| and two experimental 



groups 0, lU have now independently confirmed the 
existence of such solutions and discovered further un- 
predicted behaviour. In an annular ladder (Fig. 0a) 
and in a linear ladder |, § (Fig. 0b) breathers were 
observed with various numbers of "vertical" (or ra- 
dial) junctions rotating. The voltage drop associated 
with rotating vertical junctions causes some of the 
horizontal junctions to rotate and various configura- 
tions were commonly observed (Fig. ^). 

However the most interesting and puzzling feature 
was the switching between breathers as the bias cur- 
rent was varied. Once a breather has been initialised 
it is stable, but if the bias current is slowly decreased 
then, at some critical value of the current, the config- 
uration of rotating junctions may change thus form- 
ing a new breather. Furthermore, it appears from 
the data that the new breather may have a different 
centre of symmetry from the old breather despite the 
symmetry of the ladder itself and the symmetry of the 
driving currents about a single vertical junction. We 
present here the first numerical studies of the switch- 
ing behaviour of roto-breathers and show that this 
effect arises from the exponential growth of even in- 
finitesimal perturbations during a chaotic transient 
which is associated with the switching process. 

Aiming to describe the experiments as accurately 
as possible, we first built up a model in which the 
current / through a junction is determined by the 
well-established RSJ model equation 
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where I c is the critical current and ip — A9 — 
J A.dl where A9 is the change in superconduct- 



ing order parameter 9 across the junction and A is 
the vector potential. The "vertical" (i.e. radial) 
junctions may differ in area from the "horizontal" 
junctions by an anisotropy parameter r/ — I c h/Icv — 
C h /C v = Rv/Rh where I ch (I cv ), C h (C v ) and R h 
(Ru) are, respectively, the critical current, capaci- 
tance and resistance of a horizontal (vertical) junc- 
tion. From Fig. |]a we see that there are three un- 
knowns per plaquette: 9j, 9'j and fj. To solve for 
these unknowns we construct three equations per pla- 
quette as follows. The first two equations are ob- 
tained from current conservation at the top (inner) 
and bottom (outer) rails. The third equation is ob- 
tained by making the approximation that the induced 
flux fj — f a (where f a is the applied flux, assumed 
constant) is produced solely by the currents flowing 
around the immediate perimeter of the jth plaquette 
(Fig. 0c): 
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where 0l is an inductance parameter and we have 
chosen I cv = 1. Eq. (|2|) introduces a significant dif- 
ference from previous models jl4|, [l5| which assume 
that the induced field is proportional to the loop (or 
"mesh") current circulating the plaquette. For square 
plaquettes, Eq. (|^) represents an improvement since 
it properly takes into account the flux produced by 
the vertical conductors, although it is still not a full 
nearest- neighbour mutual inductance calculation 

Using Landau gauge we then have the following 
three coupled differential equations for each plaque- 
tte: 



_d_ 



-4sin07 + sin07 +1 + j\fj-x - ffl) 



dt 2 



dt 



sin^t 1 cosx,-_ 1 -sinx+ cosx,(4) 



9 J+1 -9J + 2r)Xj 



8tt 



{fj -fa)- 2?7sin) 



where 



x7 = 



9j — 6j 9j, 6j — 9'j + wj, A; >, , 
9J +1 - 2nfj), X j = \{0j - - 2nF/N) and N 
is the number of plaquettes. Putting F = 0, this is 
a sufficient set of equations for describing the linear 
ladder (with the imposition of the appropriate bound- 
ary conditions). For the case of the annular ladder 
we also need one further equation for the loop flux F 
(proportional to the total current circulating around 
the annulus): 
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where F a is the applied flux, assumed constant, and 
Bl is the annulus inductance parameter. Note that 
the equations remain invariant if all Oj are rotated 
by the same amount, an amount which can also be 
varied with time. It is only differences in 9 + from one 
site to another which are important. Note also that 
the equations remain invariant if 9 J and/or 9j are 
rotated by Air (not 2ir). 

In the RSJ model, on which our equations are 
based, each ideal Josephson junction has in parallel 
with it a capacitance and a resistance, or in other 
words a frequency filter. The cut-off frequency is 
1/(RC) = a. Thus voltage oscillations of frequency 
lo » a tend to be suppressed, being shorted out by 
the capacitance. In the experiments, and in our sim- 
ulations, the rotating junctions have dO/dt 3> a and 
so we can expect voltage oscillations to be suppressed 
by a large factor. We might therefore expect that all 
jObserved attractors should be characterised by hav- 
ing all rates of rotation, currents and fluxes at more 
or less fixed values. Many junctions will not rotate 
at all, in which case Eq.(|]) gives us / = 7 c sin</>. For 
those that rotate Eq.(|l|) gives us I — ad<j>/dt, i.e. the 
supercurrent and capacitative current more or less 
cancel one another. We call this the dc approx- 
imation. Using this approximation one can solve 
for the distribution of currents, voltages and fluxes 
for any configuration of rotating junctions. The ob- 
served roto-breathers are in fact quite well described 
by this approximation. 
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We now focus on determining whether or not our 
model exhibits the interesting switching and symme- 
try breaking behaviour observed in the data of Binder 
et al^, 0). Attempting to mimic the experiments, we 
use a random number generator to give 6j , and fj 
extremely small, but non-zero, initial values at t = 0. 
If instead, the initial state were taken as having all 
9j = 0'. = then we could neglect Eq.(||) and all 
solutions would maintain the symmetry xf-i = xt 
for all j; however it turns out|l4| that this symmetric 
state can be unstable to small perturbations. Let 

j . = ( Ib+Ia if J - (?) 
1 Ig otherwise 

where Ib is called the bias current. Again following 
the experiment, we slowly increase I a while keeping 
Ib = until rotation starts at site j = 0. Ia is then 
slowly decreased while at the same time increasing 
Ib to keep Iq constant. When Ib has reached the 
desired value it is then held fixed while Ia is slowly 
reduced to zero. At this point, provided Ib is not 
too large or small, we have just one rotating verti- 
cal junction (at j — 0), the horizontal junctions ro- 
tating in the /-configuration (Fig. ||). We call this 
an /-breather with Nr = 1 rotating vertical junc- 
tions. In Fig. ||a the /-breather appears to have top- 
bottom symmetry, but in fact this symmetry is bro- 
ken since 9 + (which is more or less static) actually 
varies from site to site. Also the currents circulat- 
ing the two plaquettes on either side of the rotat- 
ing junction cause a static flux pattern rather like a 
vortex-antivortex pair, the value of the flux in each 
plaquette being what one would deduce in the dc ap- 
proximation. This must cause a compressive force on 
the breather and an attractive force with the end of 
any nearby breather of the same kind. 

Finally Ib is slowly reduced while keeping I a = 0. 
The voltage current characteristics of the breathers 
produced as Ib is ramped down are shown in Fig. ||. 
We show here only the results for an annular lad- 
der with N = 8 vertical junctions and the parameter 
values otherwise chosen to mimic the annular ladder 
experiment [l?]]. We have also put Bl = since 
we find that although the inner region of the lad- 
der tends to pick up trapped flux (a vortex) in the 



breather switching process, this has little influence 
on the voltage-current characteristics for this range of 
parameter values. Note that while / and T-breathers 
occur, there are no ^-breathers. This is because in 
the annular ladder the two ends of a breather are 
connected by a quiescent region in which all super- 
conducting order parameters must be static or rotate 
together in synchrony (apart from small vibrations), 
and so Z-breathers cannot exist in the absence of 
other breathers. 

All of the breathers are well characterised by a con- 
stant resistance R (i.e. V = (d9~ /dt) = IbR) and 
the value of R agrees very closely with the value ex- 
pected in the dc approximation. Furthermore we see 
that at a critical rotation frequency the single-site 
(Nr = 1) /-breather generally broadens into a multi- 
site /-breather, but may also convert to a single-site 
or multi-site /"-breather. In agreement with the ex- 
periments, this occurs when the rotation rate of the 
horizontal junctions is close to the minimum rotation 
rate that could be supported if the junctions were 
isolated and driven by a constant current (cf. the 
miminum rotation rate of an undcrdamped pendu- 
lum driven by a constant torque). This critical fre- 
quency is more or less independent of a (for small 
a) . The multi-site /-breathers switch to T-breathers 
at the same critical frequency. In the switching, two 
of the four rotating horizontal junctions stop while 
the other two double in frequency (to maintain Kir- 
choff's Law). As Ib is reduced further the rotation 
rate of the junctions falls again until the same criti- 
cal rotation rate is reached at which point all motion 
ceases. We obtain qualitatively similar results both 
for longer (and shorter) ladders and for linear ladders 
with similar parameter values. The main difference 
in linear ladders is that Z-breathers are produced 
also. These observations agree very well with the be- 
haviour reported for a linear ladder with similar pa- 
rameter values ||. The annular ladder experiment [j] 
did not find any T-breathers but his may very well 
be due to a construction fault [fl8| . 

Ref. H shows that the linear ladder can support 
at least several types of breathers. Nearly all of the 
configurations found experimentally arose naturally 
in our simulations. Furthermore we found that in a 
multi-site breather, it appears that any two adjacent 
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vertical junctions can rotate in phase or 180° out of 
phase. It is clear that a great many breather config- 
urations are possible. 

The main result of the simulations relates to the 
switching to multi-site breathers. Although the 
single-site /-breather possesses left-right symmetry 
about j = 0, the multi-site breather produced in the 
switching process frequently has a different centre of 
symmetry. Any breather with an even number Nr 
of rotating junctions cannot have its centre of sym- 
metry at j = 0. Such breathers are produced even 
in the annular ladder (where translational invariance 
is exact) and even if we set all variables to zero ini- 
tially, in which case the dynamical equations, initial 
conditions and injected currents all have exact left- 
right symmetry about site j = 0. The origin of the 
symmetry breaking lies in the occurrence of chaotic 
dynamics during the switching. Further evidence of 
chaotic behaviour is provided by the fact that even 
the smallest change in initial conditions or of any pa- 
rameter usually results in the production of a com- 
pletely different kind of multi-site breather. Both of 
these effects (i.e. symmetry breaking and extreme 
sensitivity) were noticed in the experiments p7| . 

To study the origin of these effects we have looked 
more closely at what happens when the single-site I- 
breather becomes unstable at V — 4.3 and switching 
takes place (see Fig. ||). Approaching the switching 
instability by reducing Ib very slowly, we find that at 
1b = 0.4524 the breather appears to be stable and pe- 
riodic with no positive Lyapunov exponent. However 
when the current is lowered to Ib — 0.45235 the mo- 
tion becomes weakly chaotic with a maximum Lya- 
punov exponent (Fig. ^a) of 5 x 10~ 5 ; but the chaotic 
motion is itself unstable and switching to a new sta- 
ble breather eventually takes place after 460 000 time 
units. Thus the switching between breathers occurs 
via a chaotic transientpd^j '. Clearly we have here 
the origin of the observed symmetry breaking and 
extreme sensitivity. The nature of the resulting fi- 
nal breather is extremely sensitive to the smallest 
change in the simulation (as one might expect given 
the chaotic nature of the motion) but the time before 
switching occurs is more or less fixed. If, before wait- 
ing for switching to occur, the bias current is slowly 
lowered further to 0.45234 then the maximum Lya- 



punov exponent increases to 2 x 10 -4 and switching 
occurs after only 65 000 time units (Fig. [I]b). 

The only clear failures in the predictive power of 
the model occurred near the superconducting gap fre- 
quency where the RSJ model is bound to fail. In 
this region the experiments found that breather resis- 
tances changed and that /"-breathers could be con- 
verted to /-breathers; neither of these effects were 
reproduced by the simulations. 

We conclude that our simulations exhibit most 
of the main features of the roto-breathers recently 
observed in Josephson ladder arrays. In addition 
to reproducing the experiments we have found that 
switching occurs via a chaotic transient, and this may 
be a general feature of the switching of roto-breathers 
in Josephson arrays. This chaotic transient has some 
practical consequences. First of all it leads to the 
extreme sensitivity and symmetry breaking already 
observed in the experiments. Secondly, since it has 
has a very broad Fourier spectrum, characteristic of 
chaotic motion, we can predict the existence of a 
source of noise associated with switching. The du- 
ration of the noise (i.e. the transient) should depend 
on how quickly the bias current is ramped (Fig. ||). 
We hope that future experiments will discover the 
noise and confirm this effect. 
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Figure 1: (a) An annular ladder subjected to trans- 
verse bias currents Ij. Each circle represents a su- 
perconducting island. Each link between islands rep- 
resents a Josephson junction. 9 is the phase of the 
superconducting order parameter and / is the flux 
threading a plaquette. (b) A linear ladder, (c) Ex- 
planation of the notation used in Eq. (||). 
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Figure 2: Three configurations for the rotation of 
horizontal junctions around a breather with Nr = 1 
vertical rotating junctions: (a) the /-breather, (b) 
the T-breather, and (c) the Z-breather. 
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Figure 3: Resistance aR (where R = V/Ib) plot- 
ted against voltage V = (dO^ /dt). The results 
were obtained for an annular ladder with N = 8 
vertical junctions, a = 0.07, r\ = 0.44, /3l = 4.8, 
Bl = F a = f a = 0. The middle vertical junction was 
excited in the manner described in the text and then 
the uniform bias current Ib was ramped down. The 
data come from many runs with different starting val- 
ues of Ib- The labels T2, II etc. refer to breather 
types, e.g. "II" means the one-site /-breather. 
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Figure 4: Exponential growth of perturbations dur- 
ing the chaotic transient which occurs when the 
single-site /-breather becomes unstable at V = 
(d0-/dt) = 4.3. (a) The initial state (t = 0) was 
prepared with Ib = 0.45235. A perturbed initial 
state was then produced by introducing small ran- 
dom changes to all state variables. After an initial 
period A6q starts to grow exponentially with time. 
Switching to a new stable breather takes place at 
t = 460,000 and produces the sudden rise in A8^ 
at that time, (b) Same except the initial state was 
prepared with I B = 0.45234. 
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